Simboli i vrednosti
syms w Theta_1 Theta_2 rho_u1 R U Ug
assume(0 < R & Ug ~= 0 & 0 < U)
syms t ug
pobuda = ug == sqrt(2) * U * sin(w * t)
pobuda = 
syms R1 R2 R3 Z_C1 Z_C2
vrednosti = [R1 == 2*R/10, R2 == R, R3 == R, Z_C1 == 2*R, Z_C2 == 2*R,...
Theta_1 == pi, Theta_2 == pi/2, Ug == -1i * U]
vrednosti = 
Resavanje sistema jednacina
syms I01 I02 Id1 Id2 U01 U02 Ud1 Ud2
jednacine = [U01 == Ud1 * cos(Theta_1) + 1i * Z_C1 * Id1 * sin(Theta_1),...
I01 == 1i * 1/Z_C1 * Ud1 * sin(Theta_1) + Id1 * cos(Theta_1),...
U02 == Ud2 * cos(Theta_2) + 1i * Z_C2 * Id2 * sin(Theta_2), ...
I02 == 1i * 1/Z_C2 * Ud2 * sin(Theta_2) + Id2 * cos(Theta_2), ...
Ug == R1 * I01 + U01, Id1 == I02 + Ud1/R2, Ud1 == U02, Id2 == Ud2/R3]
jednacine = 
promenljive = [U01, I01, Ud1, Id1, U02, I02, Ud2, Id2]
promenljive = 
odziv = subs(solve(jednacine, promenljive), lhs(vrednosti), rhs(vrednosti))
odziv = struct with fields:
U01: -(U*4i)/5 I01: -(U*1i)/R Ud1: (U*4i)/5 Id1: (U*1i)/R U02: (U*4i)/5 I02: (U*1i)/(5*R) Ud2: (2*U)/5 Id2: (2*U)/(5*R)
Struja generatora
Ig = odziv.I01
Ig = 
Trenutna vrednost struje
syms ig(t)
ig(t) = simplify(subs(sqrt(2) * abs(Ig) * cos(w*t + angle(Ig)), lhs(vrednosti), rhs(vrednosti)))
ig(t) = 
Ulazna snaga
P = subs(real(Ug * conj(Ig)), lhs(vrednosti), rhs(vrednosti))
P = 
Ulazna impedansa
Zu1 = odziv.U01 / odziv.I01
Zu1 = 
Ispis ulaznog koeficijenta refleksije
subs(rho_u1 == (Zu1 - Z_C1)/(Zu1 + Z_C1), lhs(vrednosti), rhs(vrednosti))
ans =